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Summary 
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\ 

Miner's  rule  for  the  cumulative  damage  due  to  fatigue,  the  behavior 
of  which  is  well  known  in  engineering  practice  as  a  deterministic  rule, 
is  examined  from  a  probabilistic  point  of  view.  By  adopting  a  model  for 
stochastic  crack  growth  with  incremental  extensions  having  a  distribution 
with  increasing  failure  rate,  and  utilizing  some  results  from  renewal 
theory,  we  exhibit  conditions  of  dependence  upon  load  under  which  Miner's 
rule  does  yield  the  mathematical  expectation  of  the  fatigue  life.  We  also 
obtain  conditions  of  dependence  under  which  it  is  conservative  and  others 
when  it  is  unconservative.  The  relationships  between  the  mathematical 
assumptions  which  govern  when  the  rule  is,  on  the  average ,  conservative 
or  unconservative,  are  related  to  the  physical  conditions  in  practice 
which  are  known  to  force  significant  departures  from  the  rule. 


/ 


■m  *»*"  IT  II  I - 1.- 


'! 

f 

I 

-I 

i 


A 


BLANK  PAGE 

•  \ 


< 


Introduction 


In  one  of  the  first  papers  dealing  with  the  prediction  of  fatigue 
life,  M.  A.  Miner  [6]  proposed  over  21  years  ago  a  criterion  which  has 
become  famous  (or  infamous)  as  "Miner's  rule".  Dealing  with  fatigue 
as  if  it  were  a  deterministic  phenomenon  and  using  such  non-quantif iable 
concepts  as  "internal  work"  and  "damage"  the  author  derived  heuristically 
a  rule  which  has  been  applied  )  the  prediction  of  fatigue  life  under 
repeated  cyclic  loadings  with  varying  stresses. 

It  was  soon  experimentally  verified  that  Miner's  rule  could  not 
give  a  fatigue  life  prediction  under  all  types  of  loading,  which  would  be 
sufficiently  accurate  to  be  used  without  modification  as  an  engineering 
design  tool.  Thus,  since  its  formulation  in  1945  there  have  been  a  score 
or  so  other  deterministic  rules  proposed  to  take  its  place.  However, 
because  Miner's  rule  was  so  simple  to  apply  and  seemed  to  give  no  worse 
results  than  any  of  the  others,  it  has  remained  as  a  standard  for 
comparison.  In  fact,  calculations  performed  on  some  extensive  empirical 
data  [3]  showed  that  in  many  practical  cases  the  rule  yields  approximately 
correct  and  rather  conservative  results,  i.e.  that  it  tends  to  underestimate 
the  true  number  of  cycles  to  failure.  On  the  other  hand,  recent  experimental 
studies  carried  out  under  carefully  controlled  conditions,  see  [4],  [5]  and 
the  references  given  there,  have  resulted  in  data  which  make  it  seem  likely 
that  the  applicability  of  Miner's  rule  depends  on  the  arrangement  of  the 
loading  sequence  within  the  cycle,  and  that  under  some  circumstances  this 
rule  may  substantially  overestimate  the  true  number  of  cycles  to  failure. 
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2.  The  Traditional  Derivation  of  Miner's  Rule 

We  consider  standardized  specimens  of  a  material  which  are  exposed 
to  fluctuating  stresses  due  to  periodic  loading.  To  be  specific  by  a  load 
(or  load  osoillation)  we  mean  a  continuous  unimodal  function  on  the  unit 
interval,  the  value  of  which  at  any  time  gives  the  stress  imposed  by 
deflection  of  the  material  specimen.  Thus  our  load  function  contains  all 
the  information  needed  to  determine  such  parameters  as  maximum  stress, 
minimum  stress  and  average  stress  which  are  usually  used  to  define  the 
loading  oscillation. 

The  assumptions  usually  made  to  obtain  Miner's  rule  may  be 
summarized  as  follows: 

(a)  The  amount  of  damage  absorbed  by  the  material  in  any 
one  oscillation  is  determined  only  by  the  load  during 
that  oscillation. 

(b)  Each  specimen  can  absorb  the  same  amount  of  damage  and 
when  that  amount  is  attained,  failure  occurs. 

This  amount  will  be  called  the  total  damage  at  failure. 

(c)  The  total  damage  absorbed  by  the  specimen  under  a 
sequence  of  load  oscillations  is  equal  to  the  sum  of 
damages  absorbed  in  each  oscillation  during  the  sequence. 

We  introduce  the  following  notation: 

W  ■  total  damage  at  failure 

w^  ■  amount  of  damage  during  one  repetition  of  load  oscillation 

*  number  of  oscillations  to  failure  under  repeated  application 
of  the  same  load  l ^ 

where  there  are  only  i=l,...,k  possible  loads  considered. 
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Under  (a),  (b) ,  (c)  there  follows  immediately  the  equality 


(2.1) 


W  *  N^w^  f°r  each  i-l,...,k. 


If  a  sequence  of  loads  is  applied  which  are  different  from 
oscillation  to  oscillation,  and  failure  occurs  after  load  was 

applied  in  n^  oscillations  and  *n  n2»  etc.,  we  must  have 


i  niwi  =  w 
1 


and  from  (2.1)  and  (2.2)  follows 


(2.3) 


k  n. 

2  if" 1 

i-1  i 


as  a  necessary  and  sufficient  condition  for  failure  to  occur  after 
the  number  of  oscillations  to  failure  given  by 


(2.4) 


N  -  £  n  . 
i-1 


The  most  frequently  used  form  of  Miner's  Rule  is  the  following: 
If  a  given  cycle  contains  various  numbers  of  different  oscillations 
say  n^  oscillations  of  load  for  i-l,...,k  then  the  number  of 

such  cycles  that  can  be  repeated  until  failure  is 


(2.5) 


1 

k  n . 

y 

N 

i-1  i 


Hence  by  the  cumulative  damage  rule  we  determine  the  fraction  of  damage 
accrued  during  one  cycle  and  use  its  reciprocal  to  estimate  the  total 
life.  In  practice  are  determined  from  available  S-N  data  and 


JUS®*®? 


it !v*yy<»».wi»>w» . 


the  n^  are  calculated  from  a  typical  spectrum  of  loading  during  the 
cycle  and  then  we  use  (2.5)  to  estimate  N. 

The  assumptions  underlying  this  traditional  derivation  of  Miner's 
rule  are  strictly  speaking  not  verifiable  and  the  conclusion  derived, 
namely  the  deterministic  formulation  of  the  rule,  is  acceptable  only  as 
a  first  approximation.  The  following  comments  about  assumptions  (a),  (b) 
and  (c)  make  this  point  more  specific. 

Firstly  the  vague  concept  of  "damage"  is  in  need  of 

re-interpretation  in  the  light  of  modern  f ractography .  Secondly, 

Instead  of  assuming,  as  in  (a)  that  constant  repetition  of  the  same 

load  £.  on  each  oscillation  should  contribute  exactly  the  same 

amount  Wj^  of  "damage",  it  would  appear  more  plausible  to  make  the 

assumption  that  the  "damage"  caused  in  any  single  oscillation  might  vary 

from  one  oscillation  to  the  other  and  depend  upon  factors  other  than  the 

load  at  that  oscillation.  Similarly  instead  of  (b)  one  would  prefer  to 

assume  that  the  quantity  W,  the  total  "damage"  at  failure  might  also 

be  a  random  variable  which  can  assume  different  values  for  different 

specimens  of  the  material.  Finally,  assumption  (c)  that  the  "damage" 
th 

accruing  in  the  j  oscillation  is  linearly  additive  to  the  "damages" 
sustained  in  the  preceding  oscillations  appears  also  to  be  in  need  of 
re-interpretation  and  it  would  be  desirable  to  modify  it. 

In  the  next  section  these  assumptions  are  modified  in  the  sense 
suggested  above  and  some  mathematical  results  are  obtained  which  are 
analogous  to  Miner's  rule  in  its  traditional  form.  But  these  are  stated 
in  terms  of  mean  values  of  the  random  numbers  of  cycles  to  failure  and  are 
derived  under  conditions  which  appear  reasonable,  at  least  when  failure  is 
a  state  attained  during  an  early  stage  of  the  fatigue  process. 
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3.  A  Probabilistic  Model  I 

We  assume  that  fatigue  failure  is  due  to  the  growth  and  ultimate 
extension  of  a  dominant  crack.  At  each  oscillation  of  the  imposed  stress 
this  crack  is  extended  by  some  amount  which  is  a  random  function,  due  to 
the  variation  in  the  material  and  to  the  influence  of  environment,  of  the 
magnitude  of  the  imposed  stress  as  well  as  the  geometry  of  the  specimen. 
The  extension  of  the  crack  at  each  oscillation  is  therefore  a 
non-negative  random  variable  whose  distribution  may  depend  upon  several 
parameters  the  nature  of  which  we  do  not  specify  now. 

Let  be  the  sequence  of  loads  which  are  to  be  applied 

at  each  oscillation  so  that  at  the  i1"*1  oscillation  load  is  imposed. 

We  suppose  that  the  loading  is  cyclic  in  the  sense  that  for  some  m  >  1 
and  all  i-1, . . .  ,m 


and 


£jm+i 


£km+i  for  311  j  *  k< 


(3.1) 


£i+l(0) 


Ai(l). 


Hence  the  (j+l)st  cycle  is  the  loading  ^jm+i* ' ’ * ,£jm+m^ 
damage,  in  the  sense  of  extension  of  the  crack,  during  the 
is  for  j=0, 1 , 2 , . . . 


and  the  total 
cycle  of  loads 


Y...  “  X.  +•••+  X. 
j+1  jm+1  jm+m 


(3.2) 


where  is  the  microscopic  crack  extension  due  to  load 

f*  h 

in  the  iL  oscillation  of  the  (j+l)st  cycle.  We  now  set 


4^  applied 


S 

n 


Y.  +•••+  Y 
1  n 


(3.3) 


as  the  total  crack  length  at  the  end  of  n  cycles. 
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Now  we  make  the  assumption  1°  below  corresponding  to  (a)  in 

Section  2,  namely  that  loads  which  are  the  same  in  any  cycle  force  a 

crack  extension  which  has  the  same  distribution  no  matter  where  they 

occur  in  the  sequence  of  loads.  Specifically,  we  assume 

t  h 

1°  If  the  load  £  is  applied  at  the  i  oscillation,  then 
the  incremental  crack  extension  X^  is  a  random  variable 
with  a  probability  distribution  depending  only  upon  £. 

This  assumption  implies  the  independence  of  crack  extension  in  each 
oscillation  from  the  total  crack  length,  as  well  as  the  order  of  the 
loads.  While  this  may  be  realistic  in  the  early  stages  of  fatigue  crack 
growth,  it  is  not  applicable  in  many  practical  situations,  and  modified 
assumptions  will  be  proposed  in  Section  4. 

We  also  make  an  assumption  about  the  behavior  of  this  incremental 
crack  extension.  Specifically,  we  assume 

2°  For  each  1*1,2,...  the  random  variable  X^  is  non-negative 
and  has  a  distribution  with  an  increasing  failure  rate. 

We  believe  that  this  assumption  is  realistic,  at  least  during  the 
stable  phase  of  crack  growth,  by  the  following  reasoning:  a  well-accepted 
mechanism  of  crack  growth  says  that  the  longer  the  crack  propagates  during 
an  oscillation  the  more  blunt  it  becomes  (with  the  consequent  relaxation 
of  the  stress  intensity  at  the  crack  tip)  and  the  farther  into  unworked 
material  it  progresses.  Assumption  2°  is  compatible  with  such  an 
explanation  since  the  concept  of  increasing  failure  rate  applied  to 
incremental  crack  lengths  can  be  interpreted  as  meaning  that  the  longer 
the  incremental  crack  has  grown  as  a  result  of  one  oscillation,  the  more 
likely  it  will  not  propagate  a  given  distance  farther.  For  a  discussion 
of  the  definition  and  properties  of  increasing  failure  rate  (IFR) 
distributions  (or  random  variables),  see  [1]. 

Perhaps  assumption  2°  can  be  made  more  plausible  by  the  following 
argument:  Consider  molecular  bonds  holding  at  the  tip  of  the  crack  within 


{ 
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the  plastic  zone  of  the  metal.  The  stress  due  to  the  load  should 
attenuate  away  from  the  crack  tip.  It  is  reasonable  that  for  a  given 
stress  the  probability  of  the  successive  rupture  of  each  molecular  bond 
given  the  preceding  one  has  ruptured  must  decrease  since  the  stress  is 
relieved.  But  it  follows  that  this  is  equivalent  with  the  random  number 
of  bonds  broken  during  each  oscillation  being  IFR.  This  assertion  is 
made  explicit  in  Section  A  of  the  appendix. 

If  W  is  the  total  crack  length  at  which  failure  occurs,  the 
number  of  such  cycles  until  failure  is  the  integer  valued  random  variable 
N  defined  by  the  event 


[SN-1  <  S«i"' 

Note  that  failure  as  we  have  used  it  may  mean  anything  from  "catastrophic 
rupture"  to  "the  crack  is  of  such  a  length  that  it  is  inspectable" .  For 
this  and  other  reasons  one  might  consider  the  crack  length  W  which  is 
defined  as  failure  to  be  a  random  variable  with  a  given  distribution  function 
G.  Now  we  make  our  last  assumption: 

3°  The  crack  length  W  is  statistically  independent 
of  the  crack  length  for  all  n=l,2,...  . 

This  means  simply  that  knowing  the  length  w  at  which  failure  will  take 
place  has  no  influence  upon  the  behavior  of  the  crack  growth. 

The  conditional  distribution  of  N  given  W  is  for  n=l,2,... 

p  (w)  =  P [ N=n  I  W=w]  =P[S  ,<w,  S  >w] 
rn  1  n-1  n  — 

(3. A) 

p  (w)  =  H  .  (w)  -  H  (w) 
n  n-i  n 

where  H  is  the  distribution  of  S  tor  n=l,2,...  and  Hn  i  1. 
n  n  U 

Thus  the  distribution  of  N  is  given  by 


i 


i 


i 


t 
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OO  OD 

P[N<n]  =  I  'S  p  (w)dG(w)  =  f  [1-H  (w)]dG(w). 

j  ui  J  J  n 

n  J  x  n 


(3.5) 


If  we  let  v  denote  the  expected  number  of  cycles  to  failure,  then 


v  »  EN 


r  °°  r  “ 

J  1  npn(w)dG(w)  =  J  X  Hn(w)dG(w).  (3.6) 

«  n-1  „  n=0 


Notice  that  in  our  case,  from  our  assumption  of  periodic  loading,  the  Y 


for  j*l,2. 


j 

are  identically  distributed.  We  label  the  common  distribution 


F.  Then  H  is  merely  the  n-fold  convolution  of  F  with  itself.  Now  let 
n 

the  expected  crack  increment  per  cycle  be  denoted  by 


u  =  EY  «  /  xdF(x) 

J  -'o 


and  we  can  utilize  our  notation  to  state  a  well-known  result  from  renewal 

theory  which  is  fundamental  in  what  follows: 

Lemma:  If  Y^^,...  are  independently  and  identically  distributed 

non-negative  random  variables,  with  11  the  distribution  of  Y,  +•••+  Y  , 

n  In 

then  for  all  w  >  0  we  have  the  left-hand  inequality  holding 


w 

V 


1  1  I  «  (w) 
n=0 


(3.7) 


while  if  the  Yj  are  also  IFR,  then  the  right-hand  inequality  holds  as  well. 
A  complete  proof  of  this  is  given  on  page  53-54,  reference  [1]. 


If  we  denote  the  expected  value  of  W  by 


wdG(w) 


t  =  EW  = 
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and  substitute  the  results  (3.7)  into  equation  (3.6)  we  obtain  the 
fundamental  inequality  involving  the  means 

i-l<v<^.  (3*8) 

u  —  —  u 

Since  we  can  expect  in  our  application  that  is  quite  large  we 
could  interpret  (3.8)  as  equality  for  practical  purposes.  If  we  did 
so  we  would  have  what  is  indeed  Miner's  rule  in  a  different  form  as  we 
now  shew.  If  we  let  have  distribution  F(*,£^)  and 

00 

p,  =  EX  *  /  xdF(x:2,  )  i=l,...,m 

1  1  JQ  1 

be  the  expected  length  of  the  crack  extension  due  to  the  i1"^  load 
oscillation  it  then  follows  from  (3.2)  that 

m 

v  “  I  V  (3.9) 

i*l  1 

Of  course,  if  we  consider  the  derivation  above  in  the  special  case 
where  all  the  imposed  stresses  are  the  same  for  each  oscillation,  resulting 
from  the  same  load  function  S,^,  we  would  have  from  (3.8) 


T 


1  iVi 


(3.10) 


where  is  the  expected  number  of  oscillations  until  failure  under 

Substituting  (3.9),  (3.10)  into  (3.8)  and  simplifying  yields  the 
following  inequality 


1- 

1  Vi 


1  <  v  < 


—  —  m 


v  +1 
1  i 


(3.11) 


PWHMmmtmiMmmMmimmm  TiirimT'inrffi~inTr'“rrn>niifimii(>i(i..)i.,i.i  i.i... 
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In  view  of  the  assumption  of  periodicity  as  stated  in  (3.1),  there 

are  only  a  finite  number,  in  fact  not  more  than  m,  distinct  load  functions 

which  can  be  assumed  in  any  one  oscillation  of  any  cycle.  Say,  that  the 

number  of  distinct  loads  i>j  k,  then  by  virtue  of  1®  we  can  assume 

*'l*'..,*'k  are  t*nose  distinct  load  functions.  If  also  during  all  m 

oscillations  of  the  cycle  the  load  £.  occurs  n  times,  i=l,...,k, 
k 

n  *  m,  then  we  have  from  (3.11)  that 
i-1  1 


<  v  <^-J— 

k  n .  —  —  k  n , 

y  —  y  — — 

i  vi  i  v1 


(3.12) 


This  formula  is  a  direct  analog  of  (2.5)  stated  in  terms  of  mean  values  and 


in  this  form  is  precise. 


We  now  point  out  that  the  same  formula  (3.12)  holds  in  the  case  of 


randomized  stresses  with  the  proper  interpretation  of  the  notation. 

Suppose  that  during  each  oscillation  of  the  cycle  each  of  the  k  distinct 
loads  £^,...,£^  may  occur  at  random  with  probabilities  r^,...,r^, 
respectively;  then  we  obtain  by  a  similar  argument  that  (3.12)  holds  but 
n^  now  denotes  the  expected  number  of  oscillations  of  load  during 

the  cycle:  n^  =  r  jn. 

Now  we  note  that  in  both  cases,  programmed  as  veil  as  randomized  stress, 


we  can  write 


k  n. 

y  — i- 

ii  ENi 

as  an  approximate  equality  since  one  sees  that  in  fact 


(3.13) 


1.1  1 

-  -  1  <  -  <  - 

n ,  —  n,  —  _  n, 

v  _i  y  -A  y  _j_ 

^  V  ^  vi  vi+1 


? 


MRw*W,'W,*R' rt.iW> 


-11- 


Thus  under  assumptions  l8,  2°,  3°  we  see  (3.13)  is  a  correct 
re-interpretation  of  Miner's  rule  in  the  form  (2.5)  stated  in  terms 
of  expected  values.  Hence  under  any  physical  conditions  for  which 
these  assumptions  are  approximately  true  we  might  expect  good  agree¬ 
ment  on  the  average. 

It  is  necessary,  however,  to  delineate  carefully  what  has  been 
claimed  in  (3.13),  and  to  point  out  that,  if  some  assumptions  made  in 
deriving  this  formulation  of  Miner's  rule  are  not  satisfied,  (3.13) 
cannot  be  justified  and  indeed,  may  be  incorrect.  Two  such  instances 
are  now  presented  in  what  follows. 

The  usual  interpretation  of  Miner's  rule  takes  into  account  only 
the  maximum  values  of  the  loads  (in  fact,  n^  is  often  taken  as  the 
number  of  loads  with  the  same  peak  value  as  £^)  and  we  do  not  claim 
that  this  is  correct  procedure.  In  fact,  we  specifically  state  that 
it  is  the  number  of  loads  which  are  exactly  equal  to  which  are  to 

be  counted.  Thus  if  one  wishes  to  use  only  the  maximum  value  of  the 
loads  one  must  first  utilize  some  theory  to  reduce  the  actual  imposed 
load  to  its  equivalent.  For  example,  if  one  believes  that  the  expected 
incremental  crack  growth  per  oscillation  (in  some  other  terminology  this 
is  called  the  crack  growth  rate)  for  a  given  load  Z  is  proportional  to 

(max  £)a(max  Z  -  min  £)^,  (3.14) 


for  some  known  a,b  >  0,  (an  assumption  made  for  example 
the  article  by  W.  Weibull,  page  201  in  [7]),  then  we  must 
to  an  equivalent  load  Z°  which  has  minimum  stress  zero. 


in  [4]  and 
reduce  this 


> 

i 

p 

\ 

£ 

! 

I 


I 


4 

1 


4 

p's. 

*5* 


/ 


We  find 
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a  a 

max  H®  =  (max  H)a+^(max  l  -  min  2,)a+^ 

I 

and  it  is  the  loads  ( £® , . . . , 1°)  to  which  one  should  apply  Miner's 

i  m 

rule  when  taking  account  only  of  maximum  values  of  the  load. 

It  has  also  been  stated  that  all  theories,  for  which  the  expected 
crack  growth  per  oscillation  depends  only  upon  the  load  during  that 
oscillation  (such  as  in  assumption  1°)  must  fail  to  account  for  such 
empirically  demonstrated  behavior  as  crack  deceleration,  crack 
acceleration,  crack  jump  and  crack  arrest.  To  see  that  this  crticism 
is  at  least  partially  justified  we  consider  the  following  fairly 
general  model. 

We  utilize  the  evidence  given  e.g.  page  57  [4]  that  incremental 
crack  growth  takes  place  only  during  the  stress  rise  portion  of  an 
oscillation.  Suppose  we  take  as  the  expected  incremental  growth  per 
oscillation 

r1  + 

u  =  EX(£  )  =  /  hU .  (t)  ]  [  £ ' (t)  ]  dt  (3.15) 

J  J  •'Q  J  J 


,  +xifx>0  ...  ,  .  . 

where  x  =  _  , ,  and  h  is  some  non-negative  function. 

0  if  x  <  0 

It  can  be  seen  that  for  any  h  the  expected  crack  extension  during 
any  loading,  such  as  is  drawn  in  Section  B  of  the  appendix  would  be  exactly 
the  same  if  one  integrates  only  along  the  rise  portion,  in  the  manner  of 
(3.15)  whether  the  cycle  was  reversed  or  not.  This  claim  is  detailed  in 
Section  B  of  the  appendix. 

However,  experimental  evidence  such  as  reported  in  [8]  shows  that 


reversing  the  cyclic  order  alters  the  mean  crack  extension  significantly. 
This  effect  cannot  be  accounted  for  by  models  using  assumption  1°,  e.g. 
as  exemplified  by  (3.15). 
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4.  A  Probabilistic  Model  II 

It  is  clear  that  assumption  1°  of  the  preceding  section  is  the 
one  upon  which  our  results  depend  most  heavily  but  also  it  is  the  one 
that  is  the  most  open  to  question.  In  fact,  careful  experimentation 
shows  that  it  would  seem  to  be  false  under  many  conditions, 

(see  [4],  [5]).  Thus  in  order  to  obtain  a  model  which  is  consistent  with 
knowledge  gained  from  study  of  fracture  mechanics  and  metallurgy  we  make  some 
adaptations  in  1°  and  2°  but  throughout  this  section  we  assume  3°  holds. 

We  now  suppose  that  the  loading  is  cyclic  in  the  sense  of  (3.1) 
and  that 


H.  (0)  =  0  for  j  =  l,2 .  (4.1) 

jm 

and  instead  of  1°  we  propose  the  assumption 

1'  The  incremental  crack  extension  following  the 

application  of  the  i^  oscillation  is  a  random  variable 
with  a  distribution  which  depends  upon  the  size  of  the 
crack  at  the  start  of  the  cycle,  upon  the  actual 
incremental  extensions  caused  by  oscillations  which  have 
preceded  during  that  cycle,  and  a  certain  number, 

say  r,  of  the  preceding  loads. 

Note  that  the  actual  extensions  during  preceding  cycles  have  no  influence 
upon  subsequent  cycles.  This  assumption  is  based  on  the  empirical  evidence 
of  the  influence  of  the  crack  size,  the  prehistory  of  the  loading  and  the 
resulting  incremental  crack  extensions  during  those  oscillations 
(fluctuations)  of  the  cycles  (see  for  example  [3]). 


/ 
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The  loading  cycles  which  we  consider  all  have  zero  load  at  the 
start  of  each  cycle.  Thus,  at  that  time,  we  believe  there  is  no  memory 
beyond  that  of  the  actual  crack  size  and  the  most  recent  loads.  This 
is  plausible  if  each  load  fluctuation  causes  strain  energy  which  is 
dissipated  either  in  the  extension  of  the  crack  or  the  production  of 
a  plastic  zone  near  the  crack  tip  or  both.  Even  when  the  load  is 
relaxed,  there  may  remain  some  strain  hardening  or  other  effects  caused 
by  certain  loads  within  the  recent  past  which  can  still  influence  the 
crack  growth. 

Except  for  the  initial  cycles,  which  we  neglect,  the  exact  number 
determined  by  r,  we  have 

V‘Xl+"'+Xm  (4'2) 

where  each  in  the  j1"*1  cycle  has  a  distribution  with  parameters 

(Sj  ,Xj  , . . .  ,Xi_1 , X^)  where  \ ^  *  *  *  *  is  the  vector  consisting 

of  prior  r  loads.  The  indices  of  are  repeated  with  period  m  as 

indicated  in  (3.1).  As  before  S  ■  Y,  +•••+  Y  is  the  total  crack  length 

n  i  n 

at  the  end  of  the  n^  cycle. 

Actually  the  incremental  crack  extension  per  oscillation  might  well 
be  a  function  of  the  total  crack  length  up  to  that  loading  including  the 
extension  caused  by  the  previous  oscillation.  However,  our  assumption 
in  1'  is  the  more  general  since  we  consider  the  crack  length  only  up  to 
the  beginning  of  the  previous  cycle  and  an  arbitrary  dependence  upon 
the  incremental  growths  since  the  beginning  of  that  cycle  which  includes 
the  former  as  a  special  case. 
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We  now  make  the  further  assumption 


2'  For  a  given  value  of  total  crack  length  at  the  end  of 
the  cycle,  say  s ^ ,  and  the  observed  crack  extensions 
x^,...,x^  ^  in  the  cycle  prior  to  the  iC^  oscillation, 
the  incremental  crack  extension  for  the  i1"^1  oscillation 


is  X^(Sj,x^ , . . . ,x^_^ ,  \  . )  which  is  a  non-negative  IFR  random 
variable . 


This  assumption  that  is  an  IFR  random  variable  has  been  made 

before  and  its  justification  is  given  by  the  argument  in  Section  A  of  the 
appendix.  We  subsequently  consider  several  specializations  of  this  very 
general  assumption  1'  and  examine  their  consequences. 

We  first  assume 

1"  The  dependence  of  X ^  is  only  upon  A^. 

This  means  the  effect  of  the  actual  prior  incremental  crack  extensions 
during  that  cycle  and  the  total  crack  size  at  the  beginning  of  that  cycle 
have  no  influence  on  the  stochastic  behavior  of  and  can  be  neglected. 

It  follows  from  1"  that  the  incremental  extensions  are  independent  IFR 

random  variables.  This  leaves  us  with  conditions  somewhat  more  general 
than  2°. 

By  a  result  of  Barlow,  Marshall  and  Proschan  [2]  namely  that  the 
convolution  of  IFR  random  variables  is  IFR,  we  have  from  the  periodic  nature 
of  our  loading  that  are  independent  identically  distributed  IFR  random 

variables  and  thus  our  lemma  (3.7)  from  renewal  theory  applies  and  it  follows 


that 
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v ( A ) ,  the  expected  number  of  such  cycles  X  -  (£,,..., I  )  which 

1  m 

can  be  repeated  before  failure,  must  satisfy 


where  t  *  EW  and 
true  that 


— L„ 

U(X) 

P  (X)  - 


“  1  <  v(X)  <  ~yrr 
-  -  p  (X) 

EYj(X)  for  j*l,2 . 


(4.3) 


For  while  it  is 


Eyx) 


EX^X^ 


we  see  that  the  problem  of  supplying  adequate  data,  for  the  expected 
incremental  growth  during  a  load  when  that  load  has  been  preceded 
by  . . . ,  i,j_r)  becomes  magnified. 

Is  this  formula  as  far  as  we  can  go  toward  a  rederivation  of  Miner's 
rule?  As  it  stands  there  is  no  expression  possible  within  the  present 
model  which  can  enable  one  to  give  the  expected  life,  under  an  arbitrary 
periodic  load  from  the  information  contained  only  in  the  conventional 
S-N  diagram  since  the  S-N  diagram  tells  us  nothing  about  the  influence  of 
load  order.  Moreover,  unless  more  specific  information  on  the  nature  of 
the  dependence  of  r  and  X^  can  be  obtained  nothing  can  be  done  with 
(4.3). 

In  order  to  apply  Miner's  rule  and  use  the  S-N  diagram  to  calculate 
the  expected  life  the  set  A  of  all  load  sequences  of  m  elements  must 
contain  a  subset,  say  Aq,  to  which  assumptions  1°,  2°  do  apply.  Then 
the  remaining  obstacle  is  to  find  a  transformation  a  from  A  into  Aq 
such  that 


if  a(X)  -  X°,  then  p(X)  -  p(X°) 


(4.4) 


that  is,  the  random  incremental  crack  growth  per  cycle  under  A  and 
A®  both  have  the  same  expectation.  By  using  this  identity,  since 
model  I  now  applies  to  \°  in  A®,  we  have  the  same  bounds  holding 
as  in  (3.11)  and  the  same  approximation,  namely 


v(A)  - 


with  v®  *  EX(£°)  for  i=l,...,k  where  k  is  the  number  of  distinct 
loads  among  A®  and  is  the  multiplicity  of  H®  within  A®. 

Of  course  the  determination  of  the  transformation  a  is  the  real 
problem  and  the  sine  qua  non  of  a  generalized  Miner's  rule.  However, 
this  discovery  awaits  more  data  of  the  type  which  is  nuw  being 
generated  in  the  investigations  by  Schijve,  Hardrath  and  others. 

For  example  if  the  dependence  upon  X^  was  actually  a  known  function 
of  the  greatest  peak  load  in  the  preceding  r  loads  and  r  was  known  to 
be  of  moderate  value,  then  such  a  formula  could  be  developed. 


We  now  make  assumption 

1'"  The  dependence  of  X^,  the  I*"*1  incremental  extension  in 
the  cycle,  is  upon  the  prior  crack  extensions  • • • »^j_r 

and  upon  the  prior  loading  X^  *  (J^,  ^i-1’  *  '  *  ’  *i-r^  anc*  not  uPon 
the  total  crack  length  at  the  beginning  of  the  cycle. 


In  this  case  we  have  dependence  between  successive  X^'s  and 
their  convolution  need  not  be  IFR.  For  without  further  specification 
of  the  type  of  dependence,  we  can  only  conclude  that  the  are 
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independent,  non-negative  and  identically  distributed.  We  then  have 
only  half  the  fundamental  inequality  holding,  namely 


uM 


-  1  <  \j(A)  . 


Thus  if  the  conditions  are  such  that  the  dependence  of  Y^  on 
Sj_^  can  be  neglected,  perhaps  whenever  the  total  crack  length  is 
not  large,  e.g.  failure  being  defined  as  the  crack  becoming  only  large 
enough  to  be  inspectable,  then  we  see  Miner's  rule  must  always  be 
conservative.  In  case  the  transformation  a  can  be  found  and  a(\)  *  A' 
then  one  would  have 


v(A)  > 


—  k  n . 

vj« 


-  1. 


Lastly  if  1'  holds,  the  need  not  be  independent  so  that  in 

this  case  Miner's  rule  may  well  be  unconservative.  If  the  influence  of 
the  dependence  of  the  Yj  is  such  that  the  enlarging  of  the  total  crack 
size  causes  the  subsequent  random  crack  increments  per  fluctuation  to 
become  stochastically  larger  (roughly  the  probability  that  they  will  be 
larger  is  increased)  then  we  have  only  the  other  half  of  the  fundamental 
inequality  for  the  expected  number  of  cycles  to  failure,  namely 


v(A) 


This  claim  is  detailed  in  Section  C  of  the  appendix.  In  such  cases  we 
have  Miner's  rule  always  being  unconservative  in  the  sense  discussed 


here. 


In  order  to  summarize  our  conclusions  in  this  study  we  repeat 

our  basic  premise.  For  a  given  total  crack  length  at  the  end  of  the 

j^1  cycle,  say  s,  and  observed  crack  extensions  x^,...,x^  ^  in  the 

cycle  prior  to  the  i^1  oscillation,  with  the  r  preceding  load 

functions  l.  ,  the  incremental  crack  extension  for  the 

l-l  l-r 

i*”'1  oscillation  is  a  random  variable  depending  possibly  u^-on  all  these 
parameters : 


Xj ( &i » •  • • , ^i_r ,xi» • • * ■ 


(4.5) 


We  have  agreed  that  this  should  reasonably  be  assumed  to  be  an  IFR 
random  variable.  Keeping  this  in  mind  we  now  see  that  if  1°  is  true, 
all  the  arguments  of  (4.5)  may  be  omitted  except  £^.  In  this  case  the 
crack  growths  per  cycle  become  independent  identically  distributed 
IFR  random  variables  and  Miner's  rule  for  the  number  of  cycles  until 
failure  is  true  in  expectation  and  S-N  data  can  be  utilized  to  predict 
the  expected  life. 


If  1"  is  true,  namely,  that  the  arguments  of  (4.5)  may  all  be 
omitted  except  £  , ...,£^  ,  then  the  random  crack  growths  per  cycle 

are  independent  and  identically  distributed  IFR  random  variables.  In 
this  case  a  generalized  Miner's  rule  for  the  number  of  cycles  until 
failure  is  still  true  in  expectation  but  information  must  be  provided 
about  the  influence  of  load  order  for  it  to  be  correctly  utilized. 

If  1'"  is  true,  namely,  that  the  argument  s  may  be  deleted  in  (4.5), 


then  the  crack  growths  per  cycle  can  be  regarded  as  only  independent  and 
identically  distributed  random  variables  not  necessarily  IFR.  In  this  case 
Miner's  rule  for  the  number  of  cycles  until  failure  is  always  less  than  or 
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equal  to  the  true  expected  number  of  cycles  until  failure,  i.e. 
it  is  conservative. 

If  1'  is  true  and  the  dependence  upon  s  is  such  that  the  crack 
growths  per  cycle  are  dependent  and  stochastically  increasing  with  the 
total  crack  length,  then  Miner's  rule  is  unconservative  and  yields  a 
number  which  exceeds  the  true  expected  number  of  cycles  until  failure. 
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Appendix 

Section  A: 

Consider  a  macroscopic  crack  within  a  material  which,  to  fix 
ideas,  we  picture  as  follows: 


K1* 


If  ' 

£ 


For  a  given  stress  imposed  let  U  be  the  (random)  number  of  bonds  broken 

(unzippered) .  Let  q^  be  the  probability  that  the  i*"*1  bond  is  broken 

s  t 

given  that  the  (i-1)  bond  is  broken  and  hence  p^  =  1  -  is  the 

probability  the  i^  bond  is  unbroken  given  the  i-lSt  bond  is  broken.  Of 

course  we  assume  that  the  probability  that  the  i^  bond  is  broken  given  that 
st 

the  (i-1)  is  unbroken  is  zero.  Hence  in  the  multiplication  theorem  of 
conditional  probabilities 


P[U=1]  -  qxP2 

p[u-2]  -  q1q2p3 


n-1 

P[U=n-l]  =  p  TT  q,  • 

i-1 


% 

jjp 

1- 


We  now  prove  the 


Lemma:  The  are  decreasing  iff  U  is  IFR. 

Proof :  By  definition  U  is  IFR  iff  r(n)  is  increasing  where 


r(n+I) 


PfU»nl 

00 

I  P[U=k] 
k-n 


n-1 

pn  *  TT  qt 

i=l 


n-1  n 

pn  TT  1i  +  p„+l  TT  ^  +• 


n-1  n 

TJ’i  •  TTqi 


n-1  n  n 

TT  %  -  TT  qA  +  TT 


n+1 


1  *  V 


q0  +‘ 


I 


from  which  the  assertion  is  clear. | 
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Section  B 


Consider  a  loading  cycle  as  follows: 


Vvl 


f. 

f 

I 


} 

i 

i 


The  incremental  extension  due  to  the  forward  cycle,  p^,  is 


m  1 


111  HI  s*  -A. 

if  =  S  V,  =  I  /  h[  £,  (t)  ][p'  (t)  ]  dt 
j  =  l  j=l  "0  -1  J 


Integrating  only  along  the  rise  portion,  from  a^  to  ,  for  j=l,...,m, 
we  have 


b, 

!  /J 


m 


f  -  S  J  h(x)dx  =  £  [H(b  )-H(a  )] 

‘  J'1  -)  ;=!  J 


where  H  is  the  integral  of  h. 


Now  to  obtain  the  crack  extension  along  the  reversed  cycle  p^  we 
integrate  from  to  b^  j=2,...,m-l  and  a^  to  bm 


H 

J 

V 

V 

i 


I 

I 


m 


w. 


«* 


UK  Vrw  ifA 
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m-1  1 


r  =  S  /  hlVl(t)J^;  ,(t)]+dt 
j=»0  J0  m  J  m“J 


.  D 

m-l  c  m-j 

V  ■  "‘V  -  »(«!>  +  I  /  h(x)d* 

J'1  Vj+1 


m-1 


H(b  ) 
m 


-  H(3l)  +  ^  tH(bj)-H(aJ+1) ]  =  Ufl 
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The  Z^'s  are  independent  and  identically  distributed  by  F,  while  the 

Y ,'s  are  dependent  such  that  Y  depends  upon  S  =  Y,  +•••+  Y  .  The 
i  r  n+1  n  1  n 

conditional  distribution  of  Y  . .  given  S  is 

n+1  n 


P[Y  ,,<x|S  =s]  =  G(x:s). 
n+1—  n 

The  distribution  of  T  =  Z.  +•••+  Z  is  the  n-fold  convolution 

n  1  n 

of  F  with  itself.  The  distribution  of  S  is  H  where  for  n=0,l,2,... 

n  n 


(C.l) 


Hn+l(x)  ’  /„  G(K-s:s)dHn(s) 


with  Hq  =  =  1  and  we  assume  that 


(C.  2) 


C-(y:0)  =  F(y)  for  all  y  >  0 


(C.  3) 


if  s1  <  s^,  then  G(x:s^)  ^GCx^). 


I 


Section  C 


Consider  two  sequences  of  non-negative  random  variables 


^1 »Z2’  ’  ’  '  ’  ^1  ^2  ’  ’  *  * 


In  other  words  the  distribution  of  the  initial  crack  increments  per  cycle 
Z^,Y^  are  the  same  but  the  influence  of  the  total  crack  size  causes  the 
subsequent  increments  to  become  stochastically  larger  for  the  Z  process. 


We  now  prove  that  for  all  x  >  0 


I  H  (x)  <  I  F(n) (x) 


n=0 


n=0 


1: 

| 

4. 

%». 
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by  showing  that  for  all  n-1,2, 


(C.A) 


H  (x)  <  F^(x)  for  all  x  >  0, 
n  ~ 


We  do  this  by  induction.  (C.4)  is  true  for  n=l  by  (C.2).  Assume 
it  to  be  true  for  n  and  examine  the  case  for  (n+1)  .  Integrating  (C.l) 
by  parts  we  have 

X  X 

H  (x)  -  /  H  (x-y)d  G(y:x-y)  <_  /  F(n)  (x-y)d  G(n:x-y) , 

**•1  Jq  n  y  Jq  y 

by  using  (C.3)  to  show  that  G(y:x-y)  is  an  increasing  function  of  y 
and  the  induction  hypothesis  (C.4).  But  reversing  the  integration  by  parts 
we  have 

X 

H  .(x)  <  [  G(x-s:s)dFn(s)  <  Fn+1(x) 

n+1  ~  JQ 


since  G(x-s:s)  F(x-s)  by  (C.2)  and  (C.3). 
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